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Abstract. We construct the Green functions (or Feynman's propagators) for the Schrodinger equa- 
tions of the form itp t + jip xx ± tx 2 ip — in terms of Airy functions and solve the Cauchy initial 
value problem in the coordinate and momentum representations. Particular solutions of the corre- 
sponding nonlinear Schrodinger equations with variable coefficients are also found. A special case of 
the quantum parametric oscillator is studied in detail first. The Green function is explicitly given 
in terms of Airy functions and the corresponding transition amplitudes are found in terms of a 
hypergeometric function. The general case of quantum parametric oscillator is considered then in 
a similar fashion. A group theoretical meaning of the transition amplitudes and their relation with 
Bargmann's functions is established. 



1. Introduction 

In this paper we discuss explicit solutions of the Cauchy initial value problem for the one- 
dimensional Schrodinger equations 

i^ + ^±tx 2 iP = 0, r(>(x,0) = tp(x) (1.1) 

with a suitable initial data on the entire real line R. The corresponding Green functions are found 
in terms of compositions of elementary and Airy functions in the coordinate and momentum repre- 
sentations. It is well-known that the Airy equation describes motion of a quantum particle in the 
neighborhood of the turning point on the basis of the stationary, or time- independent, Schrodinger 
equation [15] . [55] . [75] . [IT] , and [49J. Here we consider an application of these functions to the 
time-dependent Schrodinger equations for certain parametric oscillator. 

It is worth noting that the Green functions for the Schrodinger equation are known explicitly 
only in a few special cases. An important example of this source is the forced harmonic oscillator 
originally considered by Richard Feynman in his path integrals approach to the nonrelativistic 
quantum mechanics [26], [27], [28], [29], and [30]; see also [43j. Since then this problem and its 
special and limiting cases were discussed by many authors; see Refs. [J3], [33], [35], [15], [19], [70] 
for the simple harmonic oscillator and Refs. [3J, [16] . [35j . [51] . [59] for the particle in a constant 
external field and references therein. 
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The case of Schrodinger equation with a general variable quadratic Hamiltonian is investigated 
in Ref. [2D]; see also [21], [22], [IS], [IE], [62J, and [BE]- Here we present a few examples that are 
integrable in terms of Airy functions. In this approach, all known exactly solvable quadratic models 
are classified in terms of solutions of a certain characterization equation. These exactly solvable cases 
may be of interest in a general treatment of the linear and nonlinear evolution equations; see [12J, 
[T7] . [18] . [19] . [39] . [12], [17] . [61] . [69] and references therein. Moreover, these explicit solutions can 
also be useful when testing numerical methods of solving the time-dependent Schrodinger equations 
with variable coefficients. The solution of the quantum parametric oscillator problem found in this 
paper is also relevant. 



2. Green Function: Increasing Case 
The fundamental solution of the time-dependent Schrodinger equation 



can be found by a familiar substitution [20 



if) = A(t) e* 5 ^'*) = 1 c i(a(t)x 2 +l3(t)xy+-y(t)y 2 ) ^ 

A/27T2/i (t) 

The real-valued functions of time a(t) , /3 (t) , 7 (t) satisfy the following system of ordinary differ- 
ential equations 

^-t + a 2 = 0, (2.3) 
f + (2.4) 

where the first equation is the special Riccati nonlinear differential equation; see, for example, [32J, 
[31], [50], [57], [58], [72] and references therein. 

The substitution 

a = — , a = — , (2.6) 

which according to Ref. [50] goes back to Jean le Round d'Alembert [23J , results in the second order 
linear equation 

fx" -tfM = 0. (2.7) 

The initial conditions for the corresponding Green function are /i(0) = and //(0) = 1/2. It is 
well-known that Eq. ( 12. 7ft can be solved in terms of Airy functions which are studied in detail; 
see, for example, [T], [3], [SO], [5E], [HI], [72] and references therein. A different definition of these 
functions that is convenient for our purposes in this paper is given in the Appendix A. 

We choose /i = (1/2) a (t) and the required Green function solution of the system is given by 

a'(t) 2 b(t) 

«0 = —7-T, Po = 7TT, 7o = — 77T, ( 2 -8) 

a (t) a (t) a (t) 
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where the Airy functions a (i) = ai (t) and b (t) = bi (t) are defined by (112. 13j) and f)12.14p . respec- 
tively. Indeed, 

dt K ' a 2 

and 



-2(a- 1 )' = 2^ = -a /5 , (2.9) 



c?7 fb\ b'a-ba' W (a, b) 1 2 



dt \a) a* a> 4* (2 ' 10) 

Thus, the Green function has the following closed form 

1 f.a'(t)x 2 -2xy + b(t)y 2 \ n 

G (x, y, t) = — ^=== exp I i — -f ^ , t>0 2.11 

V^a (t) V «W / 

in terms of elementary and Airy functions. 

It is worth noting that a more general particular solution has the form 

ip = K{x,y,t) = 1 e i{*(t)x 2 +p{t)*v+j(t)y 2 ) ; (2.12) 
y2ni[i (t) 

where \i = c\a (t) + c 2 b (t) with /x (0) = c 2 7^ 0, // (0) = c\ and 

cia 7 (t) + c 2 b' (t) Cl 

a = 77^ a (°) = — ' ( 2 - 13 ) 

cia (tj + c 2 o (t) c 2 

cia ft) + c 2 b [t) 

c 2 (3 2 (0)a(t) 

^^ 0) ~M^wUm- (2 - 15) 

This can be easily verified by a direct substitution into the system (I2.3p - fl2.5l) . 

3. Initial Value Problem: Increasing Case 

The solution of the Cauchy initial value problem 

dib 1 d 2 ib 

i -^ + 4d^ + tx ^ = ^ ^{x,0) = ip{x) (3.1) 
is given by the superposition principle in an integral form 

POO 

r(>(x,t) = / G(x,y,t) <p(y) dy, (3.2) 



where one should justify interchange differentiation and integration for a suitable initial function (p 
on R; a rigorous proof is given in Ref. 



The special case ip (y) = K (z, y, 0) of the time evolution operator ( 13. 2ft is 

/oo 
G(x,z,t) K(z,y,0) dz (3.3) 
-00 

and its inversion is given by 

POO 

G(x,y,t)=fx(0)\f3(0)\ / K{x,z,t) K*{y,z,0) dz, (3.4) 



4 



NATHAN LANFEAR AND SERGEI K. SUSLOV 



where the star denotes the complex conjugate. The familiar Euler-Gaussian-Fresnel integral 
and 



J °° e i(«*+2te) dz = ^ e -ib> i m a > 0, (3.5) 
allows to obtain the following transformation 



and its inverse 



^(t) = 2 A i(0)/i (t)(a(0) + 7o (t)), (3.6) 

a{t) = ao{t) -w£r^wy (3 ' 7) 

3(t)- f381 

/3(t) --2(a(0) +7o W)' (3 ' 8) 

7 (t) = 7 (0) - - ^ (0) - (3.9) 
W n) 4(a(0)+ 7o (*)) 

fr>w = i&wTrt bW-iW* (3 - 10) 

/i(0)/3 (0) 

a„(t) = a(t) + ,f' (t) (3.11) 
4( 7 (0) -7(t)) 

ft) - - — - — (3 12) 

Po{) ~ 2( 7 (0)-7(t))' (3 ' 12) 

T ° (i) = - tt(Q) + %& (3 ' 13) 

in the cases (13. 3ft and ( 13. 4ft . respectively. Direct calculation shows, once again, that our solutions 
( 12. 8p and (j2.13p - (!2.15p do satisfy these transformation rules. It is worth noting that the transfor- 
mation (I3.10p - (l3.13p allows to derive our Green function from any regular solution of the system 

Q-flZSD. 



4. Oscillatory Case 
A time-dependent Schrodinger equation 

l dt + Adx* tXp [ ' 

can be solved in a similar fashion by the substitution ( 12. 2ft with 

f + t + a^O, (4.2, 

^ + a/3 = 0, (4.3) 
at 

% + \f = 0. (4.4, 
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Here /i = — (1/2) a (— t) and 

«o = Kr, /3 = — — , 7o = 45 

a(—t) a(—t) a {—t) 

The Green function is 

G(x,y,t) = - r , - exp -i-^ -r^— — — ) , ^>0 (4.6) 

and the solution of the initial value problem is given by the integral (13. 2p . 

A more general particular solution has the form (I2.12p . where fi = c\a{—t) + 02b (— t) with 
fi (0) = c 2 ^ 0, fjf (0) = -ci and 

a^- C ' a, '-;': C f. a(0) = -£l, (4.7) 
cia (— tj + c 2 o (— r J c 2 

g= ; f<°> f , (4-8) 
c x a (-t) + c 2 o {—t) 

7 = 7 (Q)- ^W^-*) (49) 
7 7W 4 ( Cl a(-i) + c 2 6 (-*))■ 1 yj 

This can be easily verified by a direct substitution into the system (14 .2p - (14. 41) or with the aid of 
the transformations fl3.6|> — fl3.9[) and (I3.10I) - (I3.13I) . We leave further details to the reader. 

5. Momentum Representation 
The Schrodinger equation (12. ip takes the form 

av, 1 

l dt W 4 X ^" U [bA) 



in the momentum representation; see, for example, Ref. [22] for more details. The substitution (12.21) 
results in 

doc 1 , n „ s 

- + i -W=0, (5.2) 
^-4to/3 = 0, (5.3) 



eft 

The Riccati equation 05. 2p by the standard substitution 



dt 

- t0 2 = 0. (5.4) 



a = --^ (5.5) 
At fi ( J 

is transformed to the second order linear equation 

H" - -// -tfi = 0, (5.6) 
whose linearly independent solutions are the derivatives of Airy functions a' (t) and b' (t) . 
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We choose /i = —2b' (t) and the required solution of the system is 

b(t) a 1 a'(t) 

a °-~Wify Po ~W{ry 7 °-~46^y (5 - 7) 

The Green function is given by 

„, 1 f b (t) x 2 - 2xy + a' (t) y 2 \ , , 

G(x,y,t) = - exp -4^ — - L ^-), t>0. 5.8 

A more general particular solution has the form (12.121) . where \i = C\a! (t) + c 2 b' (t) , // (0) = c\ ^ 
and 

1 c x a (t) + c 2 b (t) ci , . 

a = -I Cl a'(t) + c 2& '(t)' tt(0) = -4^' (5 ' 9) 

c x a' (t) + c 2 6' (t) V ; 

, n , , Cl (3 2 (0)b'(t) 
7 = 7(0) + c 1 a-(t) + c 2 ^(t) - (5 - 11} 

This can be verified, once again, by a direct substitution into the system (I5.2p - (I5.4I) or with the aid 
of the transformations (EH]) (EUD and fl3J0l) -f l3l3l) . 



The oscillatory case is similar. The Schrodinger equation (14.11) in the momentum representation 
has the form 

dib d 2 ib 1 
dt dx 2 4 



i-t + t-^- -x 2 ^ = (5.12) 



and 



Here 



and 



^ + I + 4to 2 = 0, (5.13) 

dt 4 v ' 

^ + 4to/3 = 0, (5.14) 
dt 

^ + t(3 2 = 0. (5.15) 

a = 1^ (5.16) 
4t /i 

At" - -// + t/i = 0. (5.17) 

t 



The corresponding solutions are 



and 



/* = *<-*), «» = ^=|, a = -2^=5- ^° = s^f (5 ' 18) 



M = Cl a' (-£) + c 2 6' (-t) , /i (0) = c 1 ^0, (5.19) 

1 cio(-t) +ca6(-t) c 2 

a = T7 — s 777 — r, a(0) = - — , 5.20) 

4do'(-*) + c 2 6'(-t)' w 4ci' v ; 



SCHRODINGER EQUATION AND AIRY FUNCTIONS 7 

3- Cl ^°> (521) 

P ~ Cia > (_t) + C2 b> (-t) ' 

ci/3 2 (0) V (-t) 

7 = 7 , \ J -r— r- — r- . 5.22 

r r V ; cia' (-t) + c 2 6' (-t) v ; 



The Green function is given by 



. 1 h (-t) x 2 - 2xy + a' (-t) y 2 . 
G(x,y,t)= , exp [ i-± — '- — — * — —), t > 0. 5.23 



We leave further details to the reader. 



6. Gauge Transformation 



The time-dependent Schrodinger equation 

i|£= (^(P-A(x,t)f+V(x,t^, (6.1) 

where p = i~ l d/dx is the linear momentum operator, with the help of the gauge transformation 

V> = e~ i/(x 'V (6.2) 

can be transformed into a similar form 

1 % = {\ {P ~ A ' (X ' t))2 + V ' {X > t] ) * (6 ' 3) 
with the new vector and scalar potentials given by 

A- = A + ?L, V = V-%. (6.4) 

Here we consider the one-dimensional case only; see Refs. [JT] arid for more details. 
An interesting special case of the gauge transformation related to this paper is given by 

A = 0, V = -tx 2 , f = ~j (M 
A' = ~ V = -tx 2 - X -, (6.6) 



when the new Hamiltonian is 

H' = I( p -.4f + l/'=^(p+^) -ii 2 -^ (6.7) 



i(p-^) 2 + v» = I( P + ^) 2 -^^ 

7 ( P 2 + J (P x + X V) + -p J - txI ~ ~p 



Id 2 % ( d 



2x— +1) - tx 



„2 



4 dx 2 It V <9x 
and equation ( 12. ip takes the form 

.cM 1 <9 2, ?/> 9 , i ( n dib 



8 NATHAN LANFEAR AND SERGEI K. SUSLOV 

with a singular variable coefficient at the origin. Substitution (12.21) results in 



da 2 2 

— t H — a + a = 0, 

dt t 


(6.9) 




(6.10) 




(6.11) 


-, n -tfi = 0. 

ji t 


(6.12) 



where 

a 

As a result one can conclude that the time-dependent Schrodinger equation (16.81) has a solution of 
the form 

™2 , 



V> (x, t) = e-« /* / G (x, j,, t) p (y) dy, (6.13) 



where the Green function G (x, y, t) is given by (12. lip . This solution is not continuous when t — > + 
but it does satisfy the following modified initial condition 

hm e lx2/t *P (x, t) = (p (x) , (6.14) 

which reveals the structure of the singularity of the corresponding wave function at the origin. We 
leave further details to the reader. 



7. Particular Solutions of Nonlinear Schrodinger Equations 



One can find solutions of the corresponding nonlinear Schrodinger equations following Refs 
and [22] ■ For example, consider the case 

. di/j 1 d 2 ip 



+ 



dt 4 dx 2 

and look for a particular solution of the form 



+tx 2 ip = h(t) \ipr-ip, 



s > 



(7.1) 



V> = V OM) = K h {x,y,t) 



j( y a{t)x 2 +f3{t)xy+ 1 {t)y 2 +K{t)) ^ , 



constant. 



(7.2) 



Then equations (I2.3l) - (l2.5p hold with the general solution given by (I2.13p - (I2.15I) . In addition, 



~dl 



hit) 
^ (t) 



K (t) = K (0) - 



1 h(r) 
o V s '(t) 



dr. 



(7.3) 



The last integral can be explicitly evaluated in some special cases, say, when h (t) = \fi' (£) : 

A 



K (t) 



k(0) 



1 - s 



(fi 1 - 3 (t) - fi 1 -* (0)) , whens^l, 



k(0) - A In 



/i(0) 



(7.4) 



when s — 1. 



Here /x (0) 7^ 0; see [20J and [22] for more details. An example of a discontinuity of the initial data 
can be constructed by the method of Ref . [22J . Other cases are investigated in a similar fashion. 



with the Hamiltonian 
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8. Quantum Parametric Oscillator and Airy Functions 

The time-dependent Schrodinger equation for a parametric oscillator can be written in the form 

ih—=Hm (8.1) 

H _ p 2 | tw 2 (t) g ^ p=- — , (8.2) 
2m 2 % dx f 
where H is the Planck constant, m is the mass of the particle, ou (t) is the time-dependent oscillation 
frequency. The initial value problem of the form 

f)\tj ft 2 ,9 2 i]> muj 2 

m ^ = -^^ + — M + 5) X *' * (X ' 0) = $ {X) (8 ' 3) 

can be solved by the technique from the previous sections in terms of Airy functions. The substi- 
tution 

^(x,t)=e 1 ^^,r) (8.4) 



with 



results in 



mou 

r = ujt + 5, £ = ex, £= y^ft ( 8 ' 5 ) 



The Green function has the form 



(X ' y ' t} = V AmHtiiT) 6XP (*"2F ^ (T) +P (r) X2/ + 7 (r) V ^ ) ( ) 



with r = cut + 5, where // (5) = 0, // (5) = - (1/2) W (ai, bi) = 1/2 and 



fi(r) = -(ai(-5)bi(-T)-bi(-5)ai(-r)), (8.8) 

cm' (-r) 6z (-(5) - az (-6) bi' (-r) 
° lT ' at (-r) &j (-5) - ai (-5) bi (-r) ' 1 ] 

^ = 1 \h' ( • ( aw? v (8 ' 10) 

a? (— r) oz (— d) — az (— o) bi (— r) 

at 7 (-g) bi (-r) - ai (-r) bi' (-5) 

f[T) ai(-r)U(-5)-ai(-8)bi(-r)' { ' 



This can be derived with the aid of transformation fl3.10p - (l3.13l) . Thus 



1 1 2 1 

p~-ujt, a -, (3 -, 7~— (8-12) 

2 out out out 



as t —>■ + and the corresponding asymptotic formula is 

\2 
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where expression on the right-hand side is a familiar free particle propagator. The solution of the 
initial value problem (18.31) is given by 

/oo 
G(x,y,t) $(y) dy. (8.14) 
-oo 

We leave the calculation details to the reader and consider an application. 
The time-dependent quadratic potential of the form 



V(x,t) 



( l 

-muni x 2 , t < 0, 

Knu 2 (ut + 5) x 2 , < t < T, (8.15) 

-mw\ x 2 , t>T 
2 1 ' 



describes a parametric oscillator that changes its frequency from ujq to Ui during the time interval 
T. The continuity at t — and t = T defines the transition parameters uj and 5 as follows 

in terms of the initial uj$ and terminal uj\ oscillator frequencies. This model is integrable in terms 
of Airy functions with the help of the Green function found in this section as follows. 

When t < the normalized wave function for a state with the definite energy E ( n ] = huj (n + 1/2) 
is 01], H9]: 

.) „„(^„), (8 , 7) 

where H n (^ ) are the Hermite polynomials [1] , [3] , [1] , [52] , [53] , [56] , [67] , and [68] . When < t < T 
the corresponding transition wave function is given by the time evolution operator 



* n (ar, t) = U (t) *W = / G (x, y, t) ^ (y, 0) dy (8.18) 

J — oo 

with the Green function (I8.7p - fl8.lip . Finally, for t > T the wave function is a linear combination 

oo 

* n (x, t) = J2 °kn (T) ^ (x, t) (8.19) 

fc=0 

of the eigenf unctions 

corresponding to the new eigenvalues -E 1 ^ = feui (k + 1/2) with A; = 0,1,2,... . Thus function 
Cfcn (T) gives the quantum mechanical amplitude that the oscillator initially in state (uq, n) is found 
at time T in state (wi, fc) . 

For the transition period < t < T use the integral 

r 'e-*W H n (ay) dy = (A 2 - a 2 )" /2 H n ( - -^-^ V ReA 2 >0, (8.21) 
^-oo A \(A - a 2 ) / 
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which is equivalent to Eq. (30) on page 195 of Vol. 2 of Ref. [24] (the Gauss transform of Hermite 
polynomials), or Eq. (17) on page 290 of Vol. 2 of Ref. [25J. The initial wave function evolves in 
the following manner 

MM) = ?{^f\F ^ u ■ (^±^Y /2 ( 8 . 22 ) 

moo ( u p 7 \ _2 



x exp i— - a — ——z 7T- x 

' 2h \ 4 (ujI + 7 2 cj 2 ) ' 

/ mu uj 2 p 2 x 2 \ ( I ttujJq 

x exp rsMTM ) Hn W w " x 

where the time-dependent coefficients /i, a, /3, and 7 are given by equations (I8.8p - fl8.lip in terms of 
Airy functions with the argument r = out + 5 during the time interval < t < T. The asymptotics 
f l8.12p imply that \l/ n (2, t) — > (x, 0) as t — > + with the choice of principal branch of the radicals. 
A direct integration shows that 

-00 

|tf„ (x,t)\ 2 dx = 1, 0<t<T (8.23) 

by the familiar orthogonality relation of the Hermite polynomials. The normalization of the wave 
function holds also, of course, due to the unitarity of the time evolution operator. 

Then in view of the orthogonality of eigenfunctions f)8.20p the transition amplitudes are 

hfZ>{x,T)) * n (x,T) dx, (8.24) 

-DO 

where one can use another classical integral evaluated by Bailey: 

e~ x2x2 H m (ax)H n (bx) dx (8.25) 



Re A 2 > 0, 



DC 




yj (a 2 - A 2 ) ( & 2 -A 2 ) 



if m + n is even; the integral vanishes by symmetry if m + n is odd; see Refs. [10] and [44] and 
references therein for earlier works on these integrals, their special cases and extensions. 

The end result is Cfc n (T) = 0, if k + n is odd, and 



mi uj uj[3 2 . / uj 2 f3 2 ^ x x " 



u A(oj 2 + 7 2 ^ 2 ) +l \ a A(oo 2 + 7 2 ^ 2 ) 

u uj{3 2 ui . ( u 2 f3 2 -f 
X — - -s ™ ^t" hi U- 



4(w 2 + 7 2 u; 2 ) w V 4(w 2 + 7 2 cj 2 ) 
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o2 / 2 p2 \ \ -(fc+n+l)/2 

x| h , ; „ r— rr- — 2 a 



where 



4(^+ 7 2u;2) ^ 4(^+72^2) 

/ —k, -n 1 

C = (8.27) 
y (auj — 7a»i) u; 2 + (cuqu;! + «7u; 2 — f3 2 LU 2 /£j 

if + n is even. The terminating hypergeometric function is transformed as follows 

—k, —n 1 

l (1 _ fc _ n) i id + OI (8-28) 



1/2),. (1/2), ^ / -r. 



-s 



2F1 I ; -C 2 I , if k = 2r, n = 2s, 



(3/2) r (3/2) 



K 2F1 I r ' ; * ; -C a I , if fc = 2r + 1, n = 2s + 1. 



(3/2) r+s V 3/2 

It is valid in the entire complex plane; the details are given in Appendix B. The transformation 
( I8.28P completes evaluation of the Bailey integral (18.251) . 

Our function c^ n (T) gives explicitly the quantum mechanical amplitude that the oscillator ini- 
tially in state (ujq, n) is found at time T in state (u>i, k) . The unitarity of the time evolution operator 
implies the discrete orthogonality relation 

00 

J2 c *kn( T ) c kp (T) = 5 np (8.29) 

k=0 

for 2-F1 functions under consideration. The well-known orthogonal systems at this level are Jacobi, 
Kravchuk, Meixner and Meixner-Pollaczek polynomials; see, for example, [2], [3], [I], [5], [6], [7J, 
[9], [21], [37], [52], [53], [65], [66], [67] , [68], and references therein. This particular 2 F\ orthogonal 
system is reduced by the transformation (I8.28P to the Meixner polynomials. A group theoretical 
interpretation of the transition amplitudes and their relation with Bargmann's functions is discussed 
in section 10. 

In the limit T — > + , when the oscillator frequency changes instantaneously from ujq to uii, the 
transition amplitudes are essentially simplified. As a result Ck n (0) = 0, if k + n is odd, and 



c fcn (0) = i n T 



k + n + 1\ (ujquji\ 1/4 / 2 k+n+1 (ujx-ujq 



2 ) V 7r 2 7 y /c!n! (wq + cg^) V^i + ^0 / 



(fc+n)/2 



x ,fJ l"*' ; i(l + 2^^)Y (8.30) 

if + n is even. The discrete orthogonality relation (18.291) and transformation (I8.28P hold. The 
limit uj\ — > ujq is interesting from the view point of perturbation theory. 
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If the oscillator is in the ground state (uq, 0) before the start of interaction, the transition prob- 
ability of finding the oscillator in the nth excited energy eigenstate (ui,n) with the new frequency 
is given by |c 2fc+ i,o (T)\ 2 = and 



\c2k,o COT 



UJ y /UJ UJ 1 



(auJo + 7tui) 2 uj 2 + {ujquji — ceyu 2 + (3 2 u 2 /4^ 
(l/2) fe / (aujQ - ^ujif uj 2 + {ujquji + a^u 2 - (3 2 u 2 /A)' 



(8.31) 



x 



k\ 



\ {oilOq + 7CJ1) 2 uj 2 + (ujquji — ceyuj 2 + (3 2 lj 2 /4V 



where k — 0, 1, 2, ... and ^£1 l c 2fe,o CO I — 1 with the help of binomial theorem. If the oscillator 
is in the first excited state (uq, 1) , the transition probability of finding the oscillator in the nth 

|2 



excited state (a>i,n) is given by \c2k,i CO I 
|c 2 fc+i,i (T)| 2 = 



and 



3/2 



X 



(awo + 7^i) 2 ^ 2 + (^o^i — ol^iJ 2 + /3 2 u; 2 / 4) ' 
(3/2) fc / (aw - 7^i) 2 ^ 2 + (^o^i + "7^ 2 - P 2 ^ 2 /^)' 



(8.32) 



fc! 



(awo + ^oji) 2 uj 2 + (cjq^i — a'juj 2 + ^uj 2 /^)' 



where k — 0, 1, 2, ... and ^^L |c 2 fc+i,i (0| 2 = 1- These probabilities can be recognized as two special 
cases of the negative binomial distribution, or Pascal distribution, which gives the normalized weight 
function for the Meixner polynomials of a discrete variable [3], [I], [21], [H2], [53], and [68] . 

In a similar fashion, the probability that the oscillator initially in eigenstate (uo,n) is found at 
time T after the transition in state (ui, k) is given by |cfc n (T)| 2 = 0, if k + n is odd, and 



\Ckn COT 



\/ (oxuq + 7^i) 2 cj 2 + (ujquji — a7a; 2 + (3 2 u 2 /4y 



(8.33) 



X ■ 



klnlrc 



-r 2 



fc + n + 1 



x 



x 



t(auj — juji) 2 uj 2 + {ujquji + a^fu 2 — (3 2 uj 2 /4)' 
(aco>o + 7^i) 2 <^ 2 + (^0^1 — a^ou 2 + /3 2 u; 2 /4)' 



(k+n)/2 



-k, 



-n 



1F1 



:i-k-n) 



5(1+0 



if /c + n is even. The transformation (18.281) is, of course, valid but the square of the hypergeometric 
function can be simplified to a single positive sum with the help of the quadratic transformation 
( 113. ip followed by the Clausen formula (113.61) : 

\ 2 

; + (8-34) 

-n, —(k + n) /2 
- n) /2, —k — n 
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with 



(aw + l^if + (^0^1 — <y~fuj 2 + j3 2 u 2 /4:)' 

z = 

(awo — 7^1) uj 2 + (ujquji + a'jcu 2 — (3 u; 2 /4) ' 

More details are given in the next section. 



(8.35) 



Thus we have determined a complete dynamics of the quantum parametric oscillator transition 
from the initial state with the frequency ujq to the terminal one with the frequency U\ by explicitly 
solving the time-dependent Schrodinger equation with variable potential (I8.15P at all times. 



CO 



CO, 



(0 o l 











T 


t 



Figure 1 . The parametric oscillator frequency. 

9. Quantum Parametric Oscillator: General Case 

The general case of the parametric oscillator with a variable frequency of the form: 

{u , t < 0, 
w(t), 0<t<T, u = uj(0), uj{T)=u 1 (9.1) 
U U t>T 

(see Figure 1) can be investigated in a similar fashion. By the method of Ref. [2U] the corresponding 
transition Green function is given by 



G (x,y,t) 



1 



'2m fi (i) = ^ ^ x2 + f3 ^ Xy + 7 ^ V ^ 



< t < T, (9.2) 



where fi = fi (t) is a solution of the equation of motion for the classical parametric oscillator [40J, 

ft" + uj 2 (t) fi = 0, (9.3) 
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that satisfies the initial conditions a (0) =0 and // (0) = h/m. The coefficients of the quadratic 
form are 

,w m 2 \ti(t)ti'(t) J ( P >( T )) 2 /I 

The asymptotics (18.131) hold as £ — ► + ; see Ref. [63] for more details. 

A similar calculation gives the wave function (18.181) during the transition period < t < T as 
follows 

MM) = f^V /4 ; 1 (l^.Y /2 (9 . 6) 

Vvrmy ^ 2 n n\ (7 + ^0) \7 + «^o/ 
m /\ (3 2 >y "\ 2 



It satisfies the normalization condition (18.23j) . The continuity property \I/ n (x, t) — > ^/i "* (x, 0) holds 
as t —>■ + for the principal branch of the radicals. 

The transition amplitudes ( 18.241) are Ck n (T) = 0, if k + n is odd, and 



Cln(T) = r|^l)m'M r"" (2^)" /2 ( 9 .7) 

/x Ann! (7 + «u;o) m \7 + «u;o / 




(fc+n+l)/2 



where 



C = (9.8) 



yj {au - 7^i) 2 + (^o^i + 07 - /? 2 / 4 ) 2 



if /c+n is even. The transformation (18.28)) is applied and the unitarity of the time evolution operator 
implies once again the discrete orthogonality relation (I8.29j) for the 2 -Pi functions. Their relations 
with Meixner polynomials and Bargmann's functions are discussed in section 10. 

For the oscillator initially in the ground state (a>o, 0) the transition probability of finding the oscil- 
lator in the nth excited energy eigenstate (ui,n) with the new frequency is given by |c2fc+i,o (T)r = o 
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1^,0 cor 



\P\ y/uJ Q UJ\ 



(aw + + (uquji - <xy + /3 2 /4) 

(l/2) fe / (au - jujxf + (cj ^i + "7 - /5 2 / 4 )' 



(9.9) 



x 



fc! 



(au + 7Co>i) 2 + (cjo^i - «7 + /9 2 /4)' 



where = 0, 1,2, ... and ^^L |c 2 fc,o (T)\ 2 — 1- F° r the oscillator initially in the first excited state 
(uq, 1) the transition probability of finding the oscillator in the nth excited state n) is given by 
|c 2 fc,i (T)\ 2 = and 



k'2fc+i,i (t)\' 



3/2 



(aw + 7^i) 2 + (uqVi - «7 + /^ 2 /4)' 
(3/2) fc / {au - 7cux) 2 + (c^i + a 7 - /3 2 /4) : 



(9.10) 



fc! 



(aw + 7^i) 2 + (uJquJi - «7 + /3 2 /4)' 



where = 0, 1, 2, ... and Y2T=o l c 2fc,o (T)r = 1. Once again, these probabilities are two special cases 
of the negative binomial distribution, which gives the normalized weight function for the Meixner 
polynomials of a discrete variable [3], [1], [21], [52], [53], and [68J. 

In a similar fashion, the probability that the oscillator initially in eigenstate (uq, n) is found at 



time T after the transition in state (wi, fc) is given by \ck n (T)\* 

\p\ v/u^ZJT 



0, if fc + n is odd, and 



|c fcn (T)|' 



(ac^o + 7Co>i) 2 + (c^o^i - «7 + P 2 1^)' 



(9.11) 



x- 



fc + n + 1 



x 



x 



k\n\ix 

( (au - 7^i) 2 + {wqLOi + cry - /3 2 /4)' 
\ (acj + 7^x) 2 + (ujoUx - cry + /? 2 /4) 
fc, — n 
fc — n 



{k+n)/2 



2 Fl 



5 <i + <0 



if fc + n is even. The transformation (I8.28P is valid once again but the square of the hypergeometric 
function can be simplified to a single positive sum with the help of a quadratic transformation (113. II) 
followed by the Clausen formula (113.61) : 



-fc. 



-n 



iFx 



n) 



(9.12) 



-fc/2, 
'1-fc 



-n/2 
n) /2 



i + C 2 
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-k, -n, -(k + n)/2 
;/ ' 2 1 (1 -k-n)/2, -k-n' 



where 

1 + ( -2 = ( auJ o + 7^i) 2 + (up^i - «7 + /? 2 /4)^ ^ g i 

(aco>o - 7^i) 2 + {ujquji + cry - f3 2 /4) 2 
Thus substituting (I9.12p - fl9.13p into (19. lip one gets the final representation of the probability 
|cfcn {T)\ 2 in terms of a positive terminating 3 F 2 generalized hypergeometric function. 

For an arbitrary initial data in L 2 (R) : 

oo oo 

*<°> M) = £cff> *W M) , £ |c£»| a = 1, (9.14) 

n=0 n=0 

the wave function after the transition is given by 

/oo 
G{x,y,T) ^ (0) (y,0) (9.15) 
-oo 

OO /*OC 

= £d 0) / G ^y> T ) *1 0) M) 

n=0 

oo oo 

= E c i 0) ^(^ T ) = E c lM 1) (^T) 

n=0 n=0 

with 

oo 

4 1)= E C ^( T ) °n ( 9 - 16 ) 
n=0 

by (18.191) . A group theoretical interpretation will be given in the next section. The orthogonality 
property of the transition amplitudes (I8.29P implies the conservation law of the total probability 

oo „ oo 

EK1 =£K 0) f = 1 > ( 9 - 17 ) 

fc=0 n=0 

which follows, of course, from the conservation of the norm of the wave function during the transi- 
tion. 

Thus we have solved the problem of parametric oscillator in nonrelativistic quantum mechanics 
provided that the solution of the corresponding classical problem (I9.3P is known. A more convenient 
form of the transition amplitudes (19.71) will be given in the next section in terms of Bargmann's 
functions. Moreover, the quantum forced parametric oscillator can be investigated by the methods 
of Refs. GDI, H3I and Hfl. We leave the details to the reader. 



10. Group Theoretical Meaning of Transition Amplitudes 



The group theoretical properties of the harmonic oscillator wave functions are investigated in 
detail. In addition to the well-known relation with the Heisenberg-Weyl algebra of the creation and 
annihilation operators, the n-dimensional oscillator wave functions form a basis of the irreducible 
unitary representation of the Lie algebra of the noncompact group SU (1,1) corresponding to the 
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discrete positive series V 3 + ; see [51] . [IS] . [32] and [60]. In this paper we are dealing with the one- 
dimensional case only. 

Define the creation and annihilation operators: 



/ TTt ( ft d \ I TTt ( ft d \ 

° H = V 2^ [ UX ~ mdi) ' ° H = V 2^ I"" + m^J ' 
respectively, with the commutator: 

a (u) a ] {u) - a) (u) a (u) = 1 (10.2) 

and the Hamiltonian: 

H {u) = f ^ (a (u) a 1 " (u) + a) (u) a (w)) . (10.3) 
Their actions on "stationary" oscillator wave functions, 



are given by 

a^ n = s/K * ft _i, a^ n = + 1 (10.5) 

Introducing operators 

J + H = 1 (at M) 2 , J_ H = 1 (a M) 2 , (10.6) 

Jo (w) = | (a (a;) a* (w) + a* (w) a (a;)) = (u) , (10.7) 

one can easily verify the following commutation relations: 

[Jo, J±] = ±J±, [J+, J-] = -2 J . (10.8) 

For the Hermitian operators 

1 

2 v ^ ' " y 2i 

we get 

\Jxi Jy\ ^Jzi \Jyi J z\ ^Jxi \Jzi Jr] iJy (10.10) 

These commutation rules are valid for the infinitesimal operators of the non-compact group SU (1,1); 
see, for example, [UJ, [5Tj . [IS], [52], [HD], and [7T] for more details. 

One can use a different notation for the oscillator wave functions (110.41) as follows: 

f* 2n (x), if j = -3/4 and A = n + 1/4, 

jA \ ^2n+i (a;) , if j = -1/4 and A = n + 3/4, 

where n = 0, 1, 2, ... and the inequality A > j + 1 holds. The operators J± and Jo have an explicit 
form 

1 fmw 2 H 1 d \ H 
J ± = o ^ x — -z— f-T x— , ^0 = 7^ (10.12) 



J. = -(J+ + J-), J tf = -(J + -J_), J, = J , (10.9) 



2 \ ft Tko 2 dx J 2fku 
and their actions on the oscillator wave functions are given by 



J ± ^a = V(A=Fj)(A±j±l) ^,a±i, JbV-iA = A^a, (10.13) 
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whence 

J 3 4> i x = iti + i)4> j x (io.i4) 

with J 2 = Jq + Jo — J-J+ = Jq — Jo — J+J- . These relations coincide with the formulas that 
define the action of the infinitesimal operators J± and Jo of the Lie group SU (1, 1) on a basis \ j, A) 
of the irreducible representation T> 3 + belonging to the discrete positive series in an abstract Hilbert 
space (TTJ. Thus the even and odd wave functions \l/ n (x) of the one-dimensional harmonic oscillator 
form, respectively, bases for the two irreducible representations T> J + of the algebra SU (1,1) with the 
moments j = —3/4 for the even values of n and j = —1/4 for odd n. It corresponds to the double 
valued representations of the group SU (1, 1) , or quadruple valued representations of the SO (2, 1) ; 
see [52] and [60] for more details. It is important for the purpose of this paper that these group 
theoretical properties are valid "instantaneously", when u = u (t) is an arbitrary function of time. 

As a result of an elementary, but rather tedious, calculation our transition amplitudes (19.71) can 
be rewritten in the form 

Ckn (T) = T J XX , (9, r, if) (10.15) 
with the new SU (1, 1) quantum numbers 



if k = 2r, n = 2s, and 



if k — 2r + 1, n — 2s + 1. The matrix elements T 3 XX , (8, r, (p) are the so-called Bargmann functions, 
or the generalized spherical harmonics of SU (1,1) [H], [52] and [71] : 



T^(e,r,^)=e-^ v (r)e-^. (10.18) 



Here 




tiy (T) = tn /r (A +J - + l)r(A' + ,- + !) L inh r_y^ (tmh ry + V 
r(2.7 + 2)V (X- j -1)\(X' - j \ 2J \ 2) 

-A+j + 1, -A'+j + l _ 1 

2j + 2 ' sinh 2 (r/2) 

and the corresponding angles are given by 

tan* = - - ^pp^jo^m (1020) 

-27^0^! - 2^ (aj - I3 2 /A) 

tany? = — o , \ / oTT (10.21) 

(aw + 7^i) {oluq - 7a>i) + (a; ^i + «7 - P /4) (co> cji - 0:7 + p /4) 

for two Euclidean rotations, and 

tanh 2 ( T -) = ( ""°- 7 " l) 9 2+ ^ 1 + a7 - /?2/4) 2 2 (10.22) 



(au + 7CJ1) 2 + (ujqUx - 0:7 + /3 2 /4) 
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for the hyperbolic one. The branch of the radical is taken such that 

1 



sinh (r/2) 



10.23) 



The following symmetry holds: if a «-> 7 and u)q <-> Ui, then (p. It interchanges the initial and 
terminal oscillator states. 

Our formulas fjlO. 15H — fll0.23|) give a clear group theoretical interpretation of the transition am- 
plitudes for the parametric oscillator in quantum mechanics: 

CkniT) = (*«|17(D¥<P) (10.24) 

— (^ip^ g-2Wo(wi)g-j(T+ln(wi/wo))Jsg-ii(5Jo(wo)^(0) ( ^ 

= e~ lXe Z-0 7? e _i(T+ln(wi/wo))Jy V ( 5 \ e" av 
= Tl x ,(0,T,<f), 

if + n is even, in terms of the generalized spherical harmonics of the SU (1,1) algebra for the 
discrete positive series T> 3 + with j = —3/4 and j = —1/4 for the even and odd oscillator functions 
respectively. Formula (19.161) gives a transformation of the "coordinates" of a wave function from 
the old basis to the new one. 

The time evolution operator has a familiar form 

U{t)=T (9, t, <f) = e -ieM") e -i(T+M"/"o))Jy e -i<pMu ) (10.25) 

in terms of the corresponding infinitesimal operators, where 

- tJ > = j + \ x k- (10 - 26) 

The action of the hyperbolic rotation operator, or Lorentz boost, on a wave function is given by 

e~ iTj ^ (x) = e r/4 ?/> (e T/2 x) . (10.27) 

Therefore, 

t{y(r) = (V$|e^S$) (10-28) 

/oo 
(x) (e T/2 x) dx, 
-oo 

if k+n is even, which gives an integral representation for the Bargmann function under consideration; 
cf. 



Thus the SU (1,1) symmetry suggests the following algebraic form 
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up ( H 2 d mujr, n 

x exp I J — — 7 -x 

F \2Huj \2mdx 2 2 

of the time evolution operator, where the group parameters 9, r, <p are governed by the oscillator 
transition dynamics through formulas fllO.20p - fllO.23l) back to the classical equation of motion for 
the parametric oscillator fj9.3H . 

It is worth noting that Bargmann's functions are studied in detail. The functions t 3 xx , (r) are 
related to the Meixner polynomials — the unitarity property of the Bargmann functions, 

oo 

E t lx»(r)t{, x „(T) = 5 xx/ , (10.30) 

A" =7+1 

gives the discrete orthogonality relation of these polynomials. A connection with a finite set of 
Jacobi polynomials orthogonal on an infinite interval is also relevant. All basic facts about the 
functions T^ x , (6, r, (p) can be derived from the well-known properties of the Meixner and Jacobi 
polynomials; see Refs. [52], [Slj . and [7T] for more details. 

11. Summary 



dib 1 d^ib 

i 1 ^ 7 + 11 ^±tx 2 ip = (11.1) 



The time-dependent Schrodinger equations with variable coefficients 

dip 1 d 2 i 
dt 4 dx 
have the Green functions of the form 

1 / a'(±t) - 2xy + b(±t)y 2 \ 

G(x,y,t) = ^====exp(±z^^ y K )y ), t > 0, (11.2) 

^±irta{±t) \ a[±t) J 

where a (t) = ai (t) and b (t) = hi (t) are solutions of the Airy equation /i" — tfi = that satisfy the 
initial conditions a (0) = b' (0) = and a' (0) = b (0) = 1; see Appendix A below for construction of 
these solutions. 

In the momentum representation the corresponding Schrodinger equations with variable coeffi- 
cients 

dib d 2 ib 1 9 , 

have the Green functions of the form 

~, 1 ( .b(±t) -2xy + a' (±t)y 2 \ n AS 

{x - v ' t] ^ T^mm exp ( Ti wm ) • t>0 ' (1L4) 

where a' (t) = ai' (t) and b' (t) = hi' (t) are solutions of the equation fi" — ( 1 /£)// — tyt = that 
satisfy the initial conditions a' (0) = 1 and V (0) = 0; see Appendix A for further properties of these 
functions. 

Solution of the corresponding Cauchy initial value problem is given by the time evolution operator 
as follows 

/■oo 

ib(x,t)= G(x,y,t) ip(y) dy, ib (x, 0) = ip (x) (11.5) 



for a suitable function ip on R; see Ref. [S3] for more details. Additional integrable cases are given 
with the help of the gauge transformation. 
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Particular solutions of the corresponding nonlinear Schrodinger equations are obtained by the 
methods of Refs. [2D] and [22] ■ A special case of the quantum parametric oscillator with the 
Hamiltonian of the form (18.31) is studied in detail. The Green function is explicitly evaluated in 
terms of Airy functions by equations (I8.7j) — ( IHTTTT) and the corresponding transition amplitudes are 
given in terms of a hypergeometric function by formula (I8.26p . A discrete orthogonality relation for 
certain 2F1 functions is derived from the fundamentals of quantum physics. It is identified then as 
orthogonality property of special Meixner polynomials with the help of a quadratic transformation. 
An extension to the general case of parametric oscillator in quantum mechanics is also given. 
Relation of the transition amplitudes with unitary irreducible representations of the Lorentz group 
SU (1, 1) is established. Further extension to the quantum forced parametric oscillator is left to the 
reader. 

We dedicate this paper to Professor Richard Askey on his 75th birthday for his outstanding 
contributions to the area of classical analysis, special functions and their numerous applications, 
and mathematical education. 



12. Appendix A: Solutions of Airy Equation 



Bessel functions are defined as 



and the modified Bessel functions are 



_ ; o k\r(u + k + iy 

For an extensive theory of these functions, see Refs. [I], [3], [53], [56], [67], |72j and references 
therein. 

The Airy functions satisfy the second order differential equation 

u" -tu = 0. (12.3) 



Their standard definitions are 



and 



Ai(t) = ^ (I_ 1/s (z) - h /3 (*)) , (12.4) 
Bi{t) = Jl(l-i/z(z) + h/s(z)) (12.5) 



Ai(-t) = ^(j_ 1/3 (z) + J 1/3 (z)), (12.6) 
Bi(-t) = j| (J-1/3 (*) - Ji/s (*)) (12.7) 



with z = (2/3) t 3 / 2 . The Wronskian is equal to 



W(Ai{t),Bi{t)) = ~ (12.8) 
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and the derivatives are given by 

Ai'(t) = ^(l 2/3 (z)-I_ 2/3 (z)), (12.9) 
Bi'(t) = ^{l 2/3 (z)+I_ 2/3 (z)) (12.10) 

and 

Ai'(-t) = t -(j 2/3 (z)-J_ 2/3 (z)), (12.11) 
Bi'(-t) = -±=(J 2/3 (z) + J_ 2/3 (z)) (12.12) 

with z = (2/3) t 3 / 2 . 

In this paper we use the following pair of linearly independent solutions 

'2 

32/3- ^ 3 j " -v* 
~ 1 fc i (aio\ ~ 



ait) = ai(t) = ^T I ±) (12.13) 



*!(4/3) fc ^ \3j k (3k + l)\ 



k=o y ' /ft fc=0 
1 + 2 2 3 + 2 3 3 2 7 + 



and 



K*) = « (*) = (§)* 1/2 /-i/ 3 Q* 3/2 ) (12.14) 



V (t 3 /9) fc _f t m 

^Jb!(2/3), ^ ^ (3jb)I 
t 3 

+ 6 + 2 2 3 2 5 + "' 



with a (0) = 6' (0) = 0, a' (0) = b (0) = 1. Their relations with the standard Airy functions Ai (t) 
and Bi (t) are 



with the inverse 



a(t)\_l( -3V3r (1/3) 3-v«r (1/3) \ / (t) \ 
&(*) ; 2 V 3 2 / 3 r(2/3) 3^ (2/3) J \ Bi (t) J 

Ai (t) \ _ 1 / -3V«r (2/3) 3^r(l/3)Wa(t)\ 

fli(f) y< ~ tt V 3 2/3 r(2/3) 3^ (1/3) ^6(0 J l ^ ibj 



and the Wronskian is 

The derivatives are given by 



W{a(t),b{t)) = -1. (12.17) 



«'(*) = «' (t) = (0 t/_ 2/3 (V) (12.18) 



fci n /av L J U 



^*!(V3) fc V3A(3*)! 
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1+ 3" + 2332 + "' 



and 



b'(t) = W / (*) = 3^rQ)tJ 2/s Q^) (12.19) 



2 £^ ! ( 5 /3), ^ W fc (3* + 2)! 
2 + 2 2 3-5 + "' 



with the Wronskian 

W(a'(t),b'(t))=t. (12.20) 
More facts about the Airy functions can be found in Refs. p], [53], and [54] . 



13. Appendix B: Some Transformations of Hypergeometric Functions 

We derive the transformation formulas f!8.28j) as follows. In the even case k = 2r and n = 2s, use 
the quadratic transformation [3], |56j : 

/ a, b \ f 2a, 2b \ 

followed by a transformation: 

where (a) T = a (a + 1) ... (a + r — 1) = T (a + r) /T (a) , for the terminating hypergeometric func- 
tion. The reflection formula for gamma function: 

T(z)T(l-z) = (13.3) 
sin7T2; 

allows to complete the proof. 

In the odd case k = 2r + 1 and n — 2s + 1, one can use the quadratic transformation ( 113.11) for 
pure imaginary values of £ with < (Im^) 2 < 2, when the series converges. Apply the familiar 
transformation [3], [53], [56] : 

2F 1 ( a, c b ; z\ =(l-z) c - a - b 2 F 1 ( C ~°' J c °~ b ; z\ (13.4) 

to get back to a terminating hypergeometric function. The end result, namely, 

, -2r-l, -2s- 1 1 
^1 . r _ s _ 1/2 i 5 (l+*C)l (13.5) 

- r - 1/2 '- s - 1/2 ; i + c 2 

1 -r-s -1/2 ' S 
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-r, — s 



1 + C 



-r - s - 1/2 

is valid by analytic continuation in the entire complex plane. Our choice of the branch of the radical 
corresponds to the correct special value at r = s = 0. Use the transformation (113.21) and reflection 
formula f 1 1 3 . 3 j) in order to complete the proof. 

The Clausen formula [3]: 

2 



2^1 



a + b + - 
2 



3^2 



2a, 26, a + b 
a + b + ]-, 2a + 26 



and the duplication formula for gamma function: 

2 2 Z -1 



r(2z) 



■r (z\ r U + - 



7T 



;i3.6) 



;i3.7) 



have been used in sections 8 and 9. 
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